Encoding Universal Computation in the Ground States of Ising Lattices 
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We characterize the set of ground states that can be synthesized by classical 2-body Ising Hamil- 
tonians. We then construct simple Ising planar blocks that simulates efficiently a universal set of 
logic gates and connections, and hence any Boolean function. We therefore provide a new method 
of encoding universal computation in the ground states of Ising lattices, and a simpler alternative 
demonstration of the known fact that finding the ground state of a finite Ising spin glass model is 
NP complete. We relate this with our previous result about emergent properties in infinite lattices. 



I. INTRODUCTION 

The Physical Church- Turing thesis [l| provides a deep 
connection between the science of computation and the 
physical universe. It posits that the dynamics of any 
known physical system can be simulated by a Turing ma- 
chine [2| , a theoretical device that consists of a finite state 
machine together with an infinite tape. Upon reflection, 
this is a remarkable result, widely believed to be cor- 
rect. An arbitrary physical system is governed by a vast 
variety of different forces, from Coulomb interactions to 
gravity, and there is no reason, a priori, to suspect that 
all of these effects can be replicated on one particular 
machine. This presents the idea of universality: a phys- 
ical system is universal if its dynamics can be used to 
simulate any other physical system. 

The prevalence of universality in commonly studied 
systems is not only a theoretical curiosity, but also has 
consequences of practical significance. Recent results in 
computer science restrict our ability to predict the behav- 
ior of such systems. Observations of universal systems led 
to the Strong Church- Turing thesis Q, which postulates 
that a Turing Machine together with a source of random- 
ness is computationally as powerful as any other existing 
universal system. Formally speaking, we say that a task 
lies in P, or is tractable, if the task can be performed 
efficiently by a Turing Machine, i.e., the time required to 
crform it scales as a polynomial of the size of the input 
. This thesis then postulates that any problem which 
lies outside P cannot be solved with resources that scale 
polynomially with respect to the size of the problem, re- 
gardless of the method of computation used. While the 
existence of Shor's algorithm in quantum mechanics may 
provide an exception to this thesis Q, it applies to all 
current classical models of computation. 

This leads to deep insights into any universal system 
that simulates a Turing machine efficiently. Suppose such 
a system simulates a Turing machine operating on an in- 
tractable problem as input. If one could efficiently com- 
pute every physical property of this system, then one 
can use it to solve the encoded problem and therefore 
violate the Strong Church- Turing thesis. Thus, such uni- 
versal systems must necessarily exhibit properties which 
no classical algorithm can efficiently compute. 



Many other universal systems have been proposed, for 
example, logic circuits @|, the Game of Life 0], Rule 110 
8], and measurement based quantum computation [§]. 
In addition to these abstract mathematical constructs, 
many surprisingly simple physical systems capable of uni- 
versal computation have also been discovered. These in- 
clude billiard balls 10], simple dynamical systems [llj 
and the dynamics of 3-dimensional majority voting cel- 
lular automata [l2| . 

This motivates an interesting question: how simple can 
a physical system be to still exhibit universality and thus 
complex behaviour? In particular, we explore what lim- 
its can be placed on a class of Hamiltonians such that 
the evaluation of their ground states still requires the ca- 
pacity to perform universal computation. We relate this 
to the ground state decision problem: given a Hamilto- 
nian H and some number E, does there exist a state with 
energy at most El 

Interestingly, the ground state decision problem is dif- 
ficult to solve even for the simple Ising lattice, which is a 
widely used model to describe collective behaviour in di- 
verse systems, as magnetism [l3| , lattice gase s [lij , neu- 
ral activity [15| and even protein folding [lg]. While an 
efficient solution is known in the case of one dimension, F. 
Barahona showed in 1982 that the computational task is 
generally NP-complete in higher dimensions [l7|, when- 
ever some of the bonds are antiferromagnetic. Here, NP 
denotes the class of non-deterministic polynomial time 
problems; an abstract class of problems whose solutions 
can be verified, but not necessarily found, in polynomial 
time. Indeed, this connection has even allowed the engi- 
neering of spin lattice Hamiltonians whose ground states 
help model and study NP-complete problems [181 ]. 

The complexity of the ground state decision problem 
suggested that such ground states could also embed 
universal computation. Indeed, this was first proven 
with the adiabatical model of quantum computation, 
where a simple Hamiltonian with known ground state 
is adiabatically evolved to the complex Hamiltonian 
whose ground state encodes the solution to the com- 
putational problem [l9], |20| . To further simplify the 
models and make them more suitable to be recreated in 
real experiments, it has been proven that it is enough 
to consider just 2-body interactions in the Hamiltonian 
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to obtain the capability of universal computation [2l| - |24j . 

In this paper we extend those studies in the classical 
case and derive a general result on what ground state 
sets can be synthesized by a m-body Hamiltonian on a 
system of n spins. Using the circuit model of compu- 
tation, we construct simple designer circuit blocks that 
can be combined to encode a universal computer in the 
ground state of 2-body Ising Hamiltonians, in such a way 
that there is a map between any given logic circuit to 
the ground states of some Hamiltonian. This encoding, 
together with the strong Church- Turing thesis, provides 
immediate implications on the computational complex- 
ity of evaluating such ground states. Furthermore, this 
allows us to provide a simple alternative proof of Bara- 
hona's result that the ground state decision problem is 
NP-complete (l7j . 

We explore the connection of this result with the 
infinite lattice case we studied in a previous work (25[. 
We showed that there are undecidable properties in the 
infinite Ising model that give rise to emergent properties 
in the physical Ising lattice. Besides, the circuit blocks 
presented here simplify the technical parts of that work. 

This paper is organized as follows. Section [TT] intro- 
duces the required background and notation. Section ITO1 
introduced the ground synthesis problem, whilst Section 
IIVI gives an alternate proof of the universality of Ising 
ground states. Section [V] explores the consequences in 
complexity of the computational difficulty of the ground 
state problem and the the relation with the infinite case 
and emergence. Section I VII presents the main conclu- 
sions. 



II. BACKGROUND AND NOTATION 

To explore how ground states can embed universal 
computation, we first address a related practical problem 
of ground state synthesis, i.e., given a set of desired 
states, is it possible to engineer a Hamiltonian whose set 
of ground states correspond to those in the desired set? 
In particular, when reality dictates certain limits on 
the interactions available, what are the corresponding 
restrictions on the possible ground states that can be 
achieved? For example, denote the state of each spin 
by either or 1, is it possible to find a Hamiltonian 
with ground states given by {000,011,101,110}? If so, 
is it possible to engineer this Hamiltonian from Ising 
interactions? The solution to the above question gives 
us the tools to engineer a set of states that are capable 
of encoding a universal circuit. 

Let us first define the nomenclature used in this paper. 
Denote the state of each spin by either or 1. A system 
of n spins is described by a binary number b = bi . . . b n G 

, where bi G {0, 1} denotes the state of the i th spin. 
Given a state b, we make the following definitions: 



• Weight: |b| is the number of Is in b. 

• 1-Sites Ones(b) is the set of indices whose corre- 
sponding spins are 1. Ones(b) := {i : bi = 1}. 

• Descendant a is a descendent of b iff Ones(a) C 
Ones(b), i.e., the 1-sites of a are subsets of 1-sites 
of b. We write this as a partial order, a ^ b. 
Dsc(b) defines the set of all descendants of b, and 
Dsc(b, k) := {a : a ^ b, |a| = k) are all descen- 
dants of b with weight k. 

A Hamiltonian on this system is defined by a function 
H : Z2 — > M that maps each state of the system to a 
corresponding energy. A general Hamiltonian is of the 
form: 

H(b 1 ,...,b n )= (1) 

where a = a\a^ ■ ■ ■ a n G Z 2 l , cij G {0, 1}, c a are arbitrary 
constants, and the summation is taken over all binary 
strings of length n. Since we can always choose a label- 
ing of the spin states such that one of the ground state 
corresponds to 0, we assert that is a ground state of H 
(i.e. H(0) = 0) without loss of generality. 

A Hamiltonian H is m-body if it does not contain in- 
teractions involving m + 1 spins or greater, i.e: c a = 
Va such that |a| > m. The general Ising model with an 
external magnetic field is a 2-body Hamiltonian of the 
form 

// y,c,,bj,, ■ y^Mjh,. (2) 

where Cj k are the interaction energies between spins j 
and k, and Mj describes the external field at site j. 

Interaction graphs provide a convenient tool to visu- 
alize Ising Hamiltonians. Given a system of n spins, 
we associate with it a graph of n vertices where each 
spin corresponds to a single vertex. We draw an edge 
between two vertices Vi and Vj if the interaction en- 
ergy between them, Cj k is non-zero. A square Ising 
model of size N is described by an interaction graph 
with vertices Vj tk where j, k = 1, . . . , N, with edge set 
E = {(Vj >k ,Vj + i jk ), (v jtkl v jtk +i)} with j, k = 1,. . . , N+l. 

The main idea of our approach is as follows. To embed 
a binary function on two bits b ou t = /(&i,&2), we con- 
struct a Hamiltonian Hf on &i,&2, ^out with the ground 
state set 

g f = {00/(00), 01/(01), 10/(10), 11/(11)}. (3) 

We see that each element of Gf satisfy b out = /(&i,&2)- 
We define the spins in state b\ and 62 as input spins, and 
the bit in state b ou t as the output spin. We say that the 
ground state Gf encodes f . 

We can then evaluate the action of / on particular 
input, i.e., f[x,y), by introducing the external biases 
on the input spins that breaks the degeneracy of Hf 
such that the state x,yf(x,y) has lower energy than 
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the other elements of Q. For example, the Hamiltonian 
Hf{oo) = Hf + bx+b-2 would have the unique ground state 
{00/(00)}. Therefore, cooling such a system to ground 
state would allow us evaluate /(0, 0), and the computa- 
tional task of solving for a ground state of this system is 
at least as hard as evaluate /(0,0). 



III. GROUND STATE SYNTHESIS 

This motivates the problem of ground state synthesis, 
i.e., given a set of desired states, is it possible to engineer 
an m-body Hamiltonian with a coinciding set of ground 
states, and if so, how? The answer of this question can 
be directly applied to designer ground states, a set of 
ground states G / specifically designed to encode a desired 
binary function /. Should we be able to construct m- 
body Hamiltonians for arbitrary /, we can establish the 
universality of the Ising model. 

We can represent H(h) and Cb as vectors in R 2n , where 
their components are indexed by all possible values of 
b e {0, 1}". Eq. (P) implies that H(h) = Lc b , where L 
is some invertible linear map. Thus, the restriction of H 
to m-body interactions leads to a set of linear equations 
that constrain H(h). More precisely, H is an m-body 
Hamiltonian iff for each H(h) with |b| = k > m, 



ff(b)=£> p 



deDsc(b.p) 



(4) 



where a p is given by the recurrence relation (see ap- 
pendix): 



1 



1 ~~ j=l a P+j 



\b\-p 
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p = m 

1 < p < m 



(5) 



This leads immediately to constraints on the ground 
state set Q if it can be m-synthesized: 

Theorem 1 Suppose H is an m-body Hamiltonian on a 
system of n spins. For each b with |b| = k > m, define 
the sets A = {b} U Dsc(b, m - 1) U Dsc(b, m - 3) U . . . 
and B = Dsc(b, m) U Dsc(b, m — 2) U . . . Then the ground 
state set Q of H must satisfy: 



AcG^BcG 

for every b with k > in. 



(6) 



Proof: Observe that a p alternates signs for each value 
of p in Eq. ([5]), thus we can write Eq. (j4]) in the form 
£be^ c btf (b) = EbeB c b#(b). H A C Q, then the 
left hand of this equation is 0. Since H(h) > by 
assumption, it follows that the right hand side must also 
be 0, and vice versa. ■ 

This theorem immediately implies that restrictions to 
m-body Hamiltonians, for any m, will also restrict the 



sets of ground states that we can synthesize. In particu- 
lar, an m-body can only implement m-wise correlations. 
Consider for example the case of an n-body system, then 
any ground state set Q that does not satisfy 

{b : wt(b) odd} cg«{b: wt(b) even} C Q (7) 

can only be synthesized by a Hamiltonian with all n bod- 
ies interacting together. One observes that the ground 
state set corresponding to the parity function on a bi- 
nary string (i.e: /(b) = |b| mod 2) violates the above 
condition, and hence cannot be simulated by any 2-body 
Hamiltonian. Thus, we cannot simulate all binary func- 
tions directly. 

The above problem can be circumvented by introduc- 
ing ancillac, additional bits within the Ising lattice that 
are not designated as either input or output bits. For ex- 
ample, consider simulation of the NAND gate, defined by 
NAND(6i, 62) = (6i® 62)01, where all arithmetic is done 
modulo 2. Directly, a Hamiltonian H^and with ground 
state set Snand = {001, 011, 101, 110} simulates NAND. 
However, NAND can also be simulated any Hamilto- 
nian on k + 3 spins, with a ground state set of the form 
Q = {OOsool, 01s il, 10si l, llsnO}, where each Sy- de- 
note binary strings of length k. 

Now consider binary functions /, g, h, simulated by 
Hamiltonians H f, H g , Hh, with outputs 6/, b g and bh- 
The functional composition f(g(bi, 62), h(b^, 64)) on the 
four input bits bi where i = 1,...4, can be simu- 
lated by the Hamiltonian H g (b\, bi, b g ) + Hh{b^, bi, bh) + 
Hf(b g , bh, b out ), where b g and bh are introduced as ancil- 
lae. 



IV. UNIVERSALITY OF ISING GROUND 
STATES 

An arbitrary Boolean circuit that takes n input bits 
and maps them to m output bits can be decomposed 
a basic logic circuit composed of the following compo- 
nents: Wires that takes a spin as input, and copies its 
state to a neighboring spin; and NAND Gate that can 
generate all Boolean functions. These require the synthe- 
sis of the following ground state sets Gwire = {00, 11} 
and Qnand = {001,011,101,110} (In standard litera- 
ture, the FANOUT gate that copies an input bit onto 
two outputs spins is also normally required. However, 
this operation can be decomposed in spin systems into 
two wires that connect to the same input spin.) 

We the convert this to a planar circuit, that is, 
one in which no wires may intersect. This requires 
the replacement of each section where a wires inter- 
sects with a SWAP gate, SWAP(6i,6 2 ) = (b 2 ,h). We 
observe that this operation can be decomposed into 
a network of three XOR gates i.e: SWAP(6i,fe2) = 
XORi(XOR 2 (XORi(6i,6 2 ))), where XORi(6i,6 2 ) = 
(61 © 62,62) and XOR 2 (6i,6 2 ) = (61,61 ® 6 2 ). We call 
this the planar circuit representation of /. 



Therefore, we can construct a square Ising Hamilto- 
nian that synthesizes / provided there exists square Ising 
Hamiltonians that implement each of basic aforemen- 
tioned components, i.e., (1) wires, (2) NAND gates and 
(3) XOR gates. To see that each of these can be simu- 
lated by a 2-body Hamiltonian, we prove the following 
lemma: 

Lemma 2 Given a set of states Q on a system of three 
spins with 000 G Q, there exists a 2-body Hamiltonian 
that synthesizes Q if and only if 

{111} U Desc(lll, 1) C Q ^ Dcsc(lll, 2) C Q. (8) 

Proof: The forward direction is special case of Eq. ([7} 
for n — 3. To observe the converse, assume Eq. is 
true. Eq. ^ implies that H is a 2-body Hamiltonian iff 
H satisfies: 



be.4 



(9) 



where A = {111} U Desc(lll, 1) and B = Dcsc(lll,2). 
To see that Eq. ^ is true, observe that if A, B C Q 
then Eq. © is satisfied trivially. Otherwise, construct 
the Ising Hamiltonian that has assignments 



H(h) 



1 



\A/g[ 



H(d) 



1 



\b/q\' 



(10) 



for all b G A/G, d G BjQ. Here \A/G\ is the number of 
elements that lie in A but outside g. ■ 

The above lemma gives us a method to construct all 
the elements of a universal circuit from 2-body nearest 
neighbor Hamiltonians. Wires can be simulated through 
Ri = bi + b 2 - 2b 1 b 2 . Lemma [2] implies that the NAND 
can be simulated directly (to see NAND can be simu- 
lated, relabel the third bit). XOR cannot be implemented 
by the ground state of a 2-body Hamiltonian on three 
spins. However, the Hamiltonian on four spins 



HxoB,{bi,b2,bA, b Q 



: (46a -3) 
2(6i6 2 + 6 2 6 



62 + 6 ) - 46 A 
-6r6 ) + 4 (11) 



with ground states {0010,0101,1001,1100} simulates 
XOR using 6^ as an ancilla. Thus, all the above gates 
can be simulated by two-body Hamiltonians. Since these 
Hamiltonians also involve at most four spins, their inter- 
action graphs must also be planar with vertices of degree 
at most three. Thus they can all be embedded in a square 
Ising Lattice with additional ancillae (See Fig. [T]), and 
hence so can /. 

Finally, we observe that each gate can be simulated 
by a Hamiltonian on at most k spins, where k is a fixed 
number. Thus, the number of spins used to simulate / 
is at most some polynomial for the number of logic gates 
used to construct /. Therefore, the square Ising model 
can simulate an arbitrary circuit efficiently. 




(ii) 




FIG. 1. The Hamiltonian that synthesizes the XOR Gate, 
Hxor, with its corresponding interaction graph (i) can be 
embedded into a 3 x 3 square Ising Lattice, (ii) Each original 
spin bi is mapped to a set of spins Bi which are linked by 
Hwire interactions. At the ground state, all spins in each 
set Bi are of the same state, and hence behave as if they are 
a single bit. 



Theorem 3 Consider an arbitrary binary function f . 
There always exists a square Ising Hamiltonian H whose 
ground states encode f. 

The above theorem allows us to encode any logic cir- 
cuit, and thus computational task, into the ground state 
of an Ising Hamiltonian. Not only is it remarkable that 
the ground state of such simple lattices are capable of 
simulating all physical processes, but this fact also allows 
us to apply the many results of computational complex- 
ity directly onto the task for computing ground states for 
an Ising Hamiltonian. 



V. COMPUTATIONAL COMPLEXITY AND 
EMERGENCE 

Any Boolean function / can be encoded as the ground 
state of an Ising Hamiltonian Hf. Suppose now that 
/ is intractable, then the Strong Church- Turing thesis 
would necessarily imply that computing a ground state 
of Hf would also be intractable. In fact, the assertion is 
stronger. Since we can potentially encode the output of / 
in the state of any spin state, the process of determining 
the ground state of any particular spin would also be 
intractable. In this final section, we will use the above 
intuition to provide lower bounds on the computational 
difficulty of the ground state problem, i.e., finding the 
ground state of some suitable two-dimensional, nearest 
neighbor Ising Hamiltonian. 

In computational complexity Q , NP denotes the class 
of problems whose solutions can be verified, but not nec- 
essarily found, in polynomial time. It encapsulates many 
computational tasks that we would like to be able to 
solve efficiently, such as p rime factoring and the travel- 
ing salesman problem [26| . The hardest of such problems 
lie in the class NP-complete. Should any NP-complete 
problem be solved efficiently, then it could be used as a 
subroutine to efficiently solve all problems in NP and im- 
ply that P = NP. While, this remains one of the biggest 
theoretical questions in computer science, popular opin- 
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ion tends to favor that P is distinct from NP, and hence 
efficient solutions of NP-complete are unlikely. 

One particular well known NP-complete problem is the 
circuit satisfiability (CSAT) problem [27[: given a circuit 
with n input bits and a single output bit described by 
a binary function /, is there a set of inputs such that 
the output is 1? Consider a given CSAT problem with 
a circuit /. Theorem [3] implies that we can construct a 
Hamiltonian Hf together with a predefined output bit b 
such that b D = f(x) for any ground state of Hf. Since we 
can modify the any Hamiltonian by a constant without 
affecting its set of ground states, we can always choose 
Hf such that its ground state energy is 0. 

Consider the ground state decision problem, does there 
exist a state with energy at most under the Hamilto- 
nian Hf — Hf + (1 — b ou t)l The perturbation 1 — b ou t 
lifts the degeneracy in Hf such that the resulting Hamil- 
tonian H'j will have a zero energy state iff there is a set of 
inputs to / such that it outputs 1. Therefore, knowledge 
of the ground state of Hf and hence b a clearly allows 
one to solve CSAT. Therefore, the ground state decision 
problem is at least NP-hard. Furthermore, since Hf is 
a Hamiltonian on a square Ising lattice that grows at 
most polynomially with the size of the circuit, it is easy 
to check whether the energy of a given state is greater 
than 0. Thus the ground state decision problem is NP- 
complete. 

We see that the above result, originally derived by 
Barahona .17] , flows as a natural consequence of applying 
Ising lattices to solve a particular NP-complete problem. 
It is stimulating to speculate then, what other important 
results could be obtained by applying the Ising model to 
other non-trivial computational problems. The Halting 
problem [2J1 is an exciting candidate; it and its general- 
izations [281 ] prove that there exist many properties of 
Turing machines that are undecidable. Such properties 
would necessarily correspond to certain properties of the 
Ising model, and it would be interesting to see if these 
properties are physically relevant. 

Another promising avenue of research is to consider 
what the limitations on the computation of ground states 
imply about the macroscopic properties of the resulting 
Ising lattice. For example, it is easy to see how our results 
can be extended to show that computing the correlation 
length of such Ising lattices is also NP-complete. This 



leads to the concept of emergence in the infinite case in 
25], following the path established by P. Anderson in 
1972 with his celebrated paper 'More is Different' [29l j . 
where he postulated that the ground state of a spin glass 
may be non-computable. 

Emergent properties of a physical system are proper- 
ties which arise from the whole and are not deducible 
from the physical interactions of the component parts. 
In 'More Really Is Different' (25|, a special case of this 
technique was applied to show that certain macroscopic 
properties of a properly chosen, 2-dimensional, infinite 
periodic Ising lattice are emergent. That is, it is pos- 
sible to embed universal circuits within infinite periodic 
Ising lattices, such that should certain macroscopic prop- 
erties be computed, one would be able to decided whether 
a arbitrary computer program would halt. The result 
naturally motivated the question: "What would happen 
should such lattices be finite?" . In this paper we see that 
in such sceneries these emergent macroscopic properties 
are connected with the known NP-complete properties 
of finite lattice Ising spin glasses. This relation (infinite 
— > undecidable, finite —¥ NP-complete) was previously 
proved as well in planar tiling problems [30], what sug- 
gests that it could be a common feature of complex uni- 
versal systems. 



VI. CONCLUSION 

We have derived the general conditions for a desired 
set of states to be the ground state of a classical 
Hamiltonian constrained to interact with a finite number 
of spins — including 2-body interactions, i.e., the Ising 
Model. We have presented a new and simple way of 
encoding universal circuit computation in the ground 
states of Ising lattices through the construction of Ising 
blocks that implement the necessary logical gates and 
connections. This result can be immediately applied to 
derive a simple version of Barahona's original proof 17] 
that the problem of finding states on Ising Hamiltonians 
is, in general, NP-complete. 
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Appendix A: Derivation of equation (f5j) 



We first define 



3c (ai,a 2 , ■ ■ ■ ,a k ) = E a p 
P =i 



d£Dsc(6,p) 



(Al) 



Thus for any b such that ||b|j > m, we have the rela- 
tion: 

H(b) = 3c b (l,l,...,a m = 1,0,..., a* =0) 

IH-i / 

= £ £ «* 

k=m+l \d6Dsc(b,fc) 

+ <£(!,!,. ..,a m = l,0,...,a fc = 0) (A2) 



Now, we note the fact 

E #(d)=.g b (/3i,/3 2 ,...,/3 m -i,0...,0) (A3) 



d£Dsc(6,m) 



To compute j3j, Consider H(d) which has exactly m Cj 
terms of the form Cd' with ||d'|| = m. Also there exists 
" b "C m terms of the form H(d). Thus to total number 
Cd' terms is " b "C™Cj. Dividing this by the total number 
of Cd' that are descendant from b gives: 



IN 

m 



|b| 
3 



WH-j 

\h\\-m 



So that 

E H ^)=9 h c 

d£Dsc(b,m) 



k - 1 

k — m 



1 < j < m 
(A4) 



k-2 
k — m I ' ' " 



,fc-m + l,l,0...,0 



(A5) 



Substituting into Eq. (IA2 



H(b) = £ H(d)+ 

d6Dsc(6,ro) 

Qc [l — ^ 1 Cfc-mi 1 — k 2 C/s— mi • ■ • 
...,l- k - {m -VC k - m ,-{k-m),0,...,0], 

(A6) 

we eliminate the a m term in the argument of g b . By 
writing: 



E H ( d )=9, 

dGDsc(6,m-l) 

k-2 

k — (m — 1) 



fc-1 
k - (m - 1) 



,k -m + 2,l,0...,0 



(A7) 



etc, we can eliminate each of Oj, 1 < j < m recursively, 
and write out an equation for H(d) entirely from the sum 
of its descendants: 

H(b) = 3 b (ai,a 2 , . . . ,aj, . . . ,a m = 1,0, ...,0) (A8) 

with 



a m = 1 



(A9) 



k — (m — j) 
k — m 



k — (m — j) 



k - 1 

a m -u-i)( k -m- j) 



a m -2 



k — (m — j) 
k-2 



(A10) 
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substituting indices p — m — j , we get: which is the recurrence relation featured. Thus, if H 

is rra-body, then the required equation is implied. Con- 
versely, if Eq. (|A2j) is satisfied, we have: 



( k-p\ INI- 1 / \ 

a p =1 - a p+1 {k - p) - a p+2 ^ 2 )-■■■ [ Cd ) =0 Vb : ]]b|| > m (A12) 

, \ k=m+l ydSDsc(b,fe) / 

„ l<P<m (All) 

* ) which has no non-trivial solutions. 



